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Abstract 

We report on the X-ray spectral analysis of the black hole candidate XTE J1752-223 in the 2009-2010 outburst, 
utilizing data obtained with the MAXI/Gas Slit Camera (GSC), the Swift/XRT, and Suzaku, which work comple- 
mentarily. As already reported by Nakahira et al. (2010) MAXI monitored the source continuously throughout the 
entire outburst for about eight months. All the MAXI/GSC energy spectra in the high/soft state lasting for 2 months 
are well represented by a multi-color disk plus power-law model. The innermost disk temperature changed from 
~0.7 keV to -^0.4 keV and the disk flux decreased by an order of magnitude. Nevertheless, the innermost radius 
is constant at ~41 I?3.5(cos i)^ 2 km, where Z^a 5 is the source distance in units of 3.5 kpc and i the inclination. 
The multi-color disk parameters obtained with the MAXI/GSC are consistent with those with the Swift/XRT and 
Suzaku. The Suzaku data also suggests a possibility that the disk emission is slightly Comptonized, which could 
account for broad iron-K features reported previously. Assuming that the obtained innermost radius represents the 
innermost stable circular orbit for a non-rotating black hole, we estimate the mass of the black hole to be 5.5 1 ±0.28 
Mq 1)3.5(003 i)~ 2 , where the correction for the stress-free inner boundary condition and color hardening factor of 
1.7 are taken into account. If the inclination is less than 49° as suggested from the radio monitoring of transient jets 
and the soft-to-hard transition in 2010 April occurred at 1-4% of Eddignton luminosity, the fitting of the Suzaku 
spectra with a relativistic accretion-disk model derives constraints on the mass and the distance to be 3.1-55 Mq 
and 2.3-22 kpc, respectively. This confirms that the compact object in XTE J1752-223 is a black hole. 

Key words: accretion disks — black hole physics — stars: individual (XTE J1752-223) — X-rays: stars 



2 



Nakahira et al. 



[Vol. 



1. Introduction 

Galactic black hole candidates (BHCs) show various states 
characterized by their spectral shapes, temporal properties, 
and luminosities. They mostly take two major states, the 
"low/hard state" and the "high/soft state", which are referred 
to as "the hard X-ray state" and "the thermal dominant state", 
respectively, in more recent classification (see McClintock & 
Remillard 2006 and references therein) since the state of a 
BHC is not always determined by an X-ray luminosity alone. 
In the low/hard state, the X-ray energy spectra are dominated 
by a power-law component with a photon index of 1.7 and a 
high energy cutoff at ~'100 keV that shows strong short- time 
variability (Grove et al. 1998). They can be explained by ther- 
mal Comptonization of soft photons from the accretion disk by 
hot plasmas with a temperature of ~10^" K (e.g., Makishima 
etal. 2008, Takahashi et al. 2008, Shidatsu et al. 2011a). 
In a typical outburst of a transient BHC, it generally exhibits 
a spectral transition from the low/hard to the high/soft state 
through the intermediate (or very high) state when the luminos- 
ity reaches ^-^10% of the Eddingtion limit. The X-ray spectra in 
the high/soft state are characterized by a ultra-soft component 
which is considered to originate from an optically-thick and 
geometrically-thin accretion disk (so called "standard disk"; 
Shakura & Sunyaev 1973). This emission can be successfully 
described by the multi-color disk (MCD) model (Mitsuda et al. 
1984, Makishima et al. 1986) with the innermost temperature 
of ~ 1 ke V. An important signature in the high/soft state is that 
the innermost radius (?!„) of the accretion disk is kept constant, 
independent of the flux or innermost disk temperature (e.g., 
Tanaka et al. 1996). Hence, it is believed that the rjn reflects 
the innermost stable circular orbit (ISCO), which is determined 
through general relativity by the mass and angular momentum 
of the black hole; the radius of ISCO for a non-spinning black 
hole is 6i?g (i?g = GM /(? is the gravitational radius, where 
G, M, c are the gravitational constant, black hole mass, and 
light velocity, respectively). By assuming that the obtained in- 
nermost Tin to the ISCO, the black hole mass can be estimated 
from the X-ray spectrum alone for a given spin parameter 

XTE J 1752-223 was first discovered with the Proportional 
Counter Array (PCA) onboard Rossi X-ray Timing Explorer 
(RXTE) on 2009 October 23 (MJD 55127) during a Galactic 
bulge scan observation (Markwardt et al. 2009). The Gas 
Slit Camera (GSC) onboard Monitor of All-sky X-ray Image 
(MAXI; Matsuoka et al. 2009) detected the source simulta- 
neously (Nakahira et al. 2009). As described in Nakahira et 
al. (2010; hereafter Paper I), MAXI continuously monitored 
the source during the entire outburst until 2010 June. The 
MAXI light curves and hardness-intensity diagram revealed 
that the source initially stayed in the low/hard states with two 
stable flux levels for about three months, and then moved into 
the high/soft state (Paper I). A radio flare was detected si- 
multaneously with the spectral transition (Negoro et al. 2010, 
Brocksopp et al. 2010), and the proper motion of the jet ejecta 
was later found with VLBI observations (Yang et al. 2010). 
Miller- Jones et al. (2011) identified it as an approaching jet 
ejected coincidently with the X-ray state transition. From the 
observed flux of the approaching jet and its upper limit of the 
receding one, they constrained the jet speed and the inclination 



angle from the line-of-sight to be > 0.66c and < 49°, respec- 
tively, by assuming that the axes of the twin jets are aligned. 
Shaposhnikov et al. (2010) estimated the black hole mass and 
distance to be 9.6±0.9 Mq and 3.5±0.4 kpc, respectively, us- 
ing the spectral-timing correlation technique (Shaposhnikov et 
al. 2009). However, the compact object mass has not been esti- 
mated via optical mass-function technique. To firmly establish 
the nature of XTE J1752-223, however, it is quite important to 
constrain the black hole mass by an independent approach. 

In this paper, we present the results from a detailed spectral 
analysis of XTE J1752-223 utilizing the MAXI/GSC data, to- 
gether with those of Swift/XRT and Suzaku, mainly focusing 
on the spectra in the high/soft state. We then discuss the con- 
straints on the black hole mass and distance based on our re- 
sults. In the appendix, we present the current status of spectral 
calibration of the MAXI/GSC using the Crab Nebula, which is 
fully reflected in this paper and in a similar work for the black 
hole candidate GX 339-4 reported by Shidatsu et al. (201 lb). 
The spectral fitting was carried out on XSPEC version 12.6. 
Errors are quoted at statistical 90% confidence limits for a sin- 
gle parameter throughout the paper 

2. MAXI observation and analysis 

2.1. MAXI observations 

As the first astronomical mission on the International Space 
Station (ISS), MAXI has been operating since 2009 August. 
The payload carries two kinds of X-ray all-sky monitors; the 
Gas-Slit Camera (GSC; Mihara et al. 2011) and the Solid- 
state Slit Camera (SSC; Tsunemi et al. 2010, Tomida et al. 
2011). The GSC consists of Xe-gas proportional counters for 
the X-ray detector and slat collimators with a slit to localize the 
direction of the incident X-rays. The counters employ resistive- 
carbon wires for detector anodes to determine the X-ray posi- 
tion along the anode wire. Twelve identical units (refereed to 
as GSC_0, GSC_9, GSC_A and GSC_B) are assembled so 
that they instantaneously cover two rectangular field of views 
(FoVs) of 3 X 160 deg^ with an equal area. 

After the operation started in 2009 August, two counters, 
GSC_6 and GSC_9 were stopped on 2009 September 8 and 14, 
respectively, for high-voltage breakdown. Other two counters 
on GSC_A and GSC_B were stopped temporarily on September 
23, because diagnostic data indicated that they may also break 
down rather soon. Consequently, the outburst of XTE J 1752- 
223 from 2009 October to 2010 June was covered by eight GSC 
camera units out of the twelve. 

The GSC FoV is limited by slats collimators. The visibil- 
ity for a target at a given sky position changes according to 
the ISS orbital motion; hereafter we call the visibility time as 
"transit". Each transit lasts for 40-150 s, and the effective area 
of each GSC counter change due to the triangular-shaped colli- 
mator transmission function, with a peak value of 4-5 cm^. The 
MAXI/GSC scanned the direction of XTE J 1752-223 2041 
times in total from the discovery on 2009 October 23 (MJD 
55127) to 2010 June 3 (MJD 55350). The total exposure times 
effective area becomes 534.8 cm^ ksec. 
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2.2. Analysis 

For the data analysis, we used the MAXI specific analysis 
tools, which were developed by the MAXI team. We analyzed 
the GSC event data version 0.3b, which include the data taken 
by counters operated at the nominal high voltage (=1650 V) but 
excluding those of anode #1 and #2 whose energy responses 
have not been enough calibrated yet. We discarded events 
taken while the GSC FoVs were interfered by the solar panels 
and other ISS pay loads. The events detected at the anode-end 
area were also screened since the background is higher therein. 
These event were cut with a condition that the photon incident 
angle (0coi; see Mihara et al. 201 1 for the definition) is higher 
than 36°. As shown in Figure 1 we carefully extracted events 
for source and background from circular regions with a 1.5° 
radius, excluding regions within 1.7° of nearby sources, GX 
5-1, GX 9+1 and SAX J1748.9-2021 in NGC 6440. We used 
only those data when the source and background regions were 
both fully scanned in a transit. 

We performed spectral analysis of the data taken between 
MJD=55200 and MJD=55293, during which both the "hard- 
to-soft" and "soft-to-hard" transitions took place (Paper I). The 
net effective exposure was 194.5 cm^ ksec for 752 transits. We 
divided the whole dataset into groups with typical lengths of a 
few days so that the spectrum in each group has enough photon 
statistics, except for the epoch around the "hard-to-soft" state 
transition when the spectrum changed rapidly. Eventually, data 
were separated into 49 groups whose exposures times effective 
area were 0.2^-^12 cm^ ksec. The energy response matrix for 
each group was calculated by the response builder (Sugizaki et 
al. 2011). The validity of the energy response was verified 
using the Crab Nebula (see Appendix). 

2.3. Light Curves during the 2009-2010 Outburst 

Figure 2 shows the MAXI/GSC light curves of XTE J1752- 
223 during the 2009-2010 outburst in the 2-4 keV, 4-10 keV, 
and 10-20keV bands, together with the hardness ratio between 
the 4-10 keV and 2-4 keV bands. They are updated from those 
presented in Paper I, after applying the latest calibration and 
the same data screening as used for the spectral analysis (sec- 
tion 2.4 and Appendix). The MAXI/GSC first detected XTE 
J1752-223 at 15:05 (UT) on 2009 October 23 (MJD 55127) 
when the 2-20 keV X-ray intensity was 30 mCrab. It mon- 
itored the source almost uniformly except for the time when 
the FoV was close to the sun (MJD 55154-55159 and MJD 
55172-55188) and when poles of the scanning axis was close 
to the source (around MJD 55300). The light curves reveal that 
the source stayed in the initial low/hard state for the extraordi- 
nary long period of about 3 months, and took two flux-plateau 
phases meantime (Paper I). 

The "hard-to-soft" state transition occurred on MJD 55218 
and the following behavior agrees well with those of the typ- 
ical BHC outbursts (Remillard et al. 2006). The X-ray in- 
tensity was peaked at 420 mCrab on 2010 January 22 (MJD 
55218), then decayed exponentially with an e-folding time of 
34 days through the high/soft state. The "soft-to-hard" tran- 
sition started on 2010 March 30 (MJD 55285), and then the 
source returned to the low/hard state on 2010 April 6 (MJD 
55292). On 2010 June 28 (MJD 55375), the source flux fell be- 
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Fig. 1. A 2-20 keV MAXI/GSC image in the galactic coordinate 
around XTE J 1752-223. The data are accumulated from MJD 55200 
to MJD 55299. The intensity contours are shown by a log scale. The 
source and background regions for the analysis are indicated by the 
shadowed regions. 



low the MAXI/GSC detection sensitivity per day (20 mCrab). 
Thus, the total duration of the outburst became about eight 
months. 

2.4. Spectral Analysis 

For the spectral analysis, we employed the standard model 
for BHCs in the high/soft state, a multi-color disk (MCD; 
diskbb in XSPEC) model plus a power-law representing the 
hard tail. The wabs (Morrison et al. 1983) model with so- 
lar abundances by Anders et al. (1982) was applied for the 
interstellar absorption. The hydrogen column density (A^h) to- 
ward XTE J1752-223 was fixed at 0.6x10^2 cm^^^ based on 
the Swift and Suzaku results as described below. The model 
has four parameters: the innermost temperature Tin and in- 
nermost radius r■^■^-^ ' for the MCD component, with the photon 
index F and normalization at 1 keV for the power-law compo- 
nent. We first applied the model to all the data. When the MCD 
component was found to be not necessary (i.e., the 90% con- 
fidence range of the MCD normalization contains zero), only 
the power-law model was applied with its F set free. When 
the MCD component was required, we fixed F = 2.2 for all 
the spectra because it is often difficult to determine from indi- 
vidual spectrum due to poor statistics and the limited energy 
coverage of the MAXI/GSC. This photon index corresponds 
to a typical value obtained when the power-law component 
is relatively strong (>50% of the total flux in the 2-20 keV 
band). Although Shaposhnikov et al. (2010) suggest from the 



nn=\/-'Vdiskbb (i^/lOl^pc) (cos «)~ 2 , where A^diskbb is the normaliza- 
tion of the diskbb model. The distance to the source of D = 3.5 kpc and 
inclination angle i=0 are assumed. 
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in the low/hard state before MJD 55214, the intermediate state 
over MJD 55215-55218, the high/soft state over MJD 55219- 
55282, and then came back to the low/hard state after MJD 
55292 through the intermediate state between MJD 55285 and 
55292. For later discussions, we calculate the weighted aver- 
age of Tin in the high/soft state to be 41.0±2.1 km, using the 
observations on MJD 55230-55282 when the contribution of 
the power-law component was sufficiently small. 

As Steiner et al. (2009) pointed out the initial variation of 
7'in observed on MJD 55216-55218 in the intermediate state 
may be caused by the ignoring of Comptonized photons in the 
calculation of Tin. In these epochs, the fraction of the power- 
law component was >50% of the total X-ray flux in the 2-20 
keV band, and hence the Compton up-scattering of the disk 
blackbody emission can significantly reduce the normaliza- 
tion of the direct MCD component (McClintock & Remillard 
2006). We find it difficult, however, to obtain a reliable esti- 
mate of the innermost radius corrected for this effect from the 
MAXI/GSC data alone, which strongly couples with the pho- 
ton index of the Compton scattered component. We leave de- 
tailed investigation of the spectra in the intermediate state for 
future work, and concentrate on those in the high/soft state in 
the following analysis. 

3. Swift Observations and Analysis 



Fig. 2. MAXI/GSC and Swift/XRT light curves of XTE J1752-223 in 
ttie 2009-2010 outburst. Four panels from the top show the GSC light 
curves in three energy bands of 2-4, 4-10 and 10-20 keV, and the hard- 
ness ratio between the 4—10 keV and 2-A keV bands. The bottom shows 
the Swift/XRT light curve in the 0.3-10 keV band. 

RXTE/PCA data that the power-law slope is variable between 
r « 1 .6 — 2.2 by using a spectral model different from ours, we 
confirmed that varying T value within this range only changes 
the best-fit MCD normalization by < 3%, and hence does not 
affect our conclusions. 

The MCD plus power-law model or single power-law model 
gave good fits for all the MAXI/GSC data. Table 1 summarizes 
the obtained spectral parameters, while figure 3 plots the evolu- 
tion of the model parameters. Figure 4 shows typical response- 
unfolded i/Fjy spectra, together with their best-fit models. 

The spectra during MJD 55200-55214 (DatalD = 1-10) re- 
quired no MCD component. The photon indices and fluxes was 
almost constant meanwhile at «1.7 and 1.2xl0~^ erg cm^^ 
s-\ respectively. On MJD 55215.64-55216.15 (DatalD = 11- 
12), the spectrum dramatically softened to F ^2.0. The MCD 
component then appeared after MJD 55216.65 (DatalD = 13) 
and lasted until MJD 55286 (DatalD = 46), although the power- 
law flux was dominant MJD 55216.65-55218.37 and on MJD 
55286.97 in the 2-20 keV band. The MCD flux reached to 
a peak of --1 x lO^^erg cm'^ on 2010 January 23 (MJD 
55219), and then decreased. The innermost temperature Tin 
gradually decreased from -^0.7 keV to ^^0.4 keV toward MJD 
55286.97. By contrast, the innermost radius was almost 
constant at -^41 km (for D = 3.5 kpc and i = 0°) then, ex- 
cept for the epoch before the state transition on MJD 55218.37 
(DatalD = 13-15). 

Based on these results, we identify that the source was likely 



From 2009 October 25 to 2010 July 29, the Swift/XRT 
(Burrows et al. 2005) carried out 67 pointing observations 
of XTE J 1752-223 in the Windowed Timing (WT) mode. Due 
to the Sun angle constraint, the observations were interrupted 
between 2009 November 4 (MJD 55139) and 2010 February 3 
(MJD 55230). We analyzed data of 17 observations taken in 
the high/soft state after 2010 February 4 (MJD 55231). All the 
XRT spectra and light curves were produced by the web inter- 
face (Evans et al. 2009) supplied by the UK Swift Science 
Data Centre at the University of Leicester. We used the XRT 
response matrix file for the WT mode version 12. All the 17 
observations with an exposure time of a few ksec each were 
analyzed independently. 

First results of the Swift/XRT spectral analysis were already 
reported by Curran et al. (2010). They employed the same 
MCD plus power-law model as we used in the MAXI/GSC 
analyses, but derived a different result that the innermost ra- 
dius changed significantly with time. The difference is consid- 
ered to come from the fact that they left the power-law index 
free and obtained a wide range of F between ^0 and ~3, even 
though this quantity is not constrained when the MCD compo- 
nent is strong. Hence we reanalyzed the XRT spectra by fixing 
F at 2.2, the same value as adopted for the MAXI/GSC spec- 
tra. The A'h was left free because the Swift/XRT is sensitive 
down to ^0.3 keV. Resultant parameters are summarized in ta- 
ble 2, which are also plotted on figure 3. The values of were 
thus constant, with at most 5% scatter around weighted mean 
of r,n=41A ±0.3 km. 
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Fig. 3. Evolution of the spectral parameters of XTE J1752-223. The MAXI/GSC, Swift/XRT, and Suzaku results using the MCD plus power-law model 
are shown by black, red, and blue points, respectively. From top to bottom panels, 2-20 keV fluxes of MCD and power-law components, the hydrogen 
column density A'^jj; the photon index T, the inner most temperature Tin, and the innermost radius rjn are indicated. 



4. Suzaku Observations and Results 

4.1. Observation 

We also observed XTE J1752-223 in the high/soft state with 
Suzaku (Mitsuda et al. 2007) as a Target-of-Opportunity (ToO) 
from 2010 February 24 (MJD 55251) 04:58:00 to February 25 
04:27:24 (ObsID: 904005010). Suzaku carries 4 X-ray tele- 
scopes (XRT, Serlemitsos et al. 2007), each with a focal-plane 
X-ray CCD imager called XIS (X-ray Imaging Spectrometer; 
Koyama et al. 2007) covering the 0.2-12 keV energy band. 
Since XIS2 has not been available, the two remaining front- 
illuminated (FI) CCDs (XISO and XIS3) were operated for 
this observation together with the back-illuminated (BI) CCD 
(XISl). The three XIS cameras were operated with burst and 
window options. The burst time, window size, and editing 
mode were "0.3 sec, 1/4 window, and 2x2" for XISO and 
XIS3, and "0.1 sec, normal window, and 3x3" for XISl, re- 
spectively. The Hard X-ray Detector (HXD; Takahashi et al. 
2007; Kokubun et al. 2007), covering the 10-70 keV energy 



band with Si PIN photo-diodes (HXD-PIN), and the 50-600 
keV range with GSO scintillation counters (HXD-GSO), was 
operated in the normal mode. The source was observed at the 
HXD nominal position. 

4.2. Data Reduction 

The data reduction and analyses were performed using 
Suzaku FTOOLS in HEASOFT version 6.10 provided by the 
NAS A/GSFC Suzaku GOF All the XIS and HXD data were re- 
processed by the Suzaku pipeline processing version 2.4. 12.27, 
with the latest calibration database (CALDBVER) hxd20090902, 
xis20090925, and xrt20080709. 

The image degradation due to "thermal wobbling" of the 
satellite pointing was corrected by using the attitude file up- 
dated with AEATTCOR2. We screened the XIS data under the 
standard criteria, and discarded events during time intervals 
when the telemetry of the XIS was saturated. Resultant net 
exposure was 10.6 ksec and 449 sec for FIs and BI, respec- 
tively. The BI exposure was short, because 0.1 sec burst time 
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Fig. 4. Examples of the i/Fu spectrum of XTE J 1752-223 observed with the MAXI/GSC together with the best-lit model. Panels correspond to DatalDs 
6, 10, 11, 14, 16, 26, 35, 44, and 49. Data ID 35 was taken approximately at the same time as the Suzaku observation. 



with normal window make the exposure 0.1 sec in 8 sec pe- 
riod and the exposure is one twelfth of each FI's exposure. 
Furthermore, the fraction of out-of-time events (Koyama et al. 
2007) is more significant in the BI data, which would intro- 
duce systematic uncertainties in the spectral analysis. Hence, 
we did not use the XISl data in this paper The average 0.5- 
10 keV count rate was ^400 cts s^^ and ^700 cts s^^ with 
FIs and BI, respectively. We extracted the XISO and 3 events 
from a rectangle region of 8f6 x4f5 tracing the 1/4 window area 
centered on the image peak, and the background from another 
region avoiding contaminating point sources. The source was 
so bright that the XIS suffered from photon pileup at the image 
center. Using the software of AEPILEUPCHECKUP (Yamada et al. 
201 1 in preparation), which automatically detects the extent of 
pileup utilizing radial surface brightness profiles and other sets 
of information, we excluded a circular within 1 '.5 of the image 
peak for XISO and 3, to extract the events whose pileup frac- 
tion is below 1%. We combined the spectra and response files 
of the two Fl cameras to increase the photon statistics in the 
spectra. To account for possible uncertainties in the instrument 
calibration (e.g. Makishima et al. 2008), systematic errors of 
1 % was assigned to each energy bin of the XIS spectra. 



We obtained the cleaned HXD events with the standard cri- 
teria. After the dead-time correction, the net exposures of PIN 
and GSO were 21.2 and 27.0 ksec, respectively. The dead- 
time fraction was 18.2 and 7.5 % for PIN and GSO, respec- 
tively. The cosmic X-ray background (CXB) was modeled 
assuming an exponentially cutoff power-law model (Boldt 
1987). The non-X-ray background (NXB) model is provided 
by the HXD team (Fukazawa et al. 2009). We used the 
model with METHOD="LCFlTDT(bgdjd)", and the version of 
METHODV="2.0ver0804" and "2.4ver09 12-64" for the data of 
PIN and GSO, respectively. After subtracting the modeled 
NXB and CXB, the source signals were significantly detected 
by PIN up to 50 keV, above the systematic uncertainty on the 
NXB model ('-^3%). Since the GSO signals were comparable 
to the systematic errors of the NXB model there in (^1%), we 
consider the source detection with GSO insignificant. 

4. 3. Spectral Analysis 

As shown in figure 5(a), we used an energy range of the 1.2- 
10 keV for the XIS spectrum and 15-50 keV for the HXD-PIN 
spectrum, where the calibrations of the energy responses were 
well established. The energy bands of 1.6-2.0 keV and 2.2-2.4 
keV were excluded to avoid large systematic uncertainties of 
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the effective area near the Silicon K edge (1.74 keV and 1.84 
for Ka and K/3, respectively) and gold M edge (2.29 keV), re- 
spectively. In the combined fit to the XIS and HXD-PIN spec- 
tra, we employed a cross normalization factor of 1.18 between 
XIS and PIN (SUZAKU-MEMO-2008-06). 

We firstly fitted the spectra with the MCD plus power-law 
model modified by interstellar absorption (wabs model), the 
same one as adopted to fit the MAXI/GSC and Swift/XRT 
spectra in the previous sections. The residuals in units of x 
are shown in the figure 5(b). The fit was found to be un- 
acceptable ix^/dof = 230.1/136), leaving an emission-line 
feature at «6.5 keV, and significant broad edge-like residu- 
als above ^1 keV. These features may be explained away by 
several different approaches; e.g., by invoking a broad Fe-K 
line concept (Reis et al. 2010) or applying so-called smeared- 
edge model (e.g., Ebisawa et al. 1994; Dotani et al. 1998) 
to account for the edge structure. The latter approach with 
wabs*smedge*(diskbb+powerlaw+gaus) model gave an ac- 
ceptable fit (x^/dof = 141.6/133; fig 5(c)) with a maximum 
optical depth of S.SOlg yg for fixed edge energy and width of 
7.11 keV and 10 keV, respectively. A narrow emission line 
was found at i?c=6.59^o 07 keV with a la line width of 10 
eV (fixed) and an equivalent width of EW=63t.ti These 
approaches are empirical and may be degenerate, however. 

As an alternative possibility, we resorted to the idea of 
Kolehmainen et al. (2011) that the disk emission is broader 
than the simplest MCD model, which could lead to artifi- 
cial detection of broad iron-K features. Assuming that weak 
Comptonization of the MCD emission took place, we replaced 
the diskbb model with a compPS model (Poutanen & Svensson 
1996), which computes Compton scattering by hot corona us- 
ing exact numerical solution radiative transfer equation. We 
chose "slab" (geom=l) geometry and a MCD spectrum as the 
seed photons. We left the optical depth of the corona r as a free 
parameter but fixed the electron temperature at To = 10 keV 
and covering fraction at covfj-ac = 1 to avoid strong coupling 
between these parameters; the choice of To did not affect our 
results on the disk parameters over the statistical errors. The 
presence of a fluorescence iron-K line at Ri6.6 keV indicates 
that a reflection continuum must be included in the model as 
well. Hence, we adopted pexnv model (Magdziarz & Zdziarski 
1995) for the power-law component, which was considered to 
be the main source of irradiation onto the disk in the high/soft 
state. The model thus became wabs*(compPS+pexnv+gaus) 
in the XSPEC terminology. For the pexriv component, we as- 
sumed no high energy cutoff, and fixed the solid angle of the 
reflector at = 27r, inclination at 25°, and ionization param- 
eter at = 1000 erg cm s^^, which consistently explain both 
iron-K line intensity and center energy. The fit was found to 
be acceptable with /dof = 129.9/132; as shown in the fig- 
ure 5(d), the broad edge-like residuals mostly disappeared. The 
model gave an absorbed 2-20 keV flux of (3.76±0.02)x 10^^ 
erg s-^cm-^. We obtained t = 0.41 ±0.05 and T = 2.13±0.01 
in the pexriv component, a reasonable value in the high/soft 
state (e.g., Zdziarski et al. 2004). More detailed analysis, more 
advanced disk emission modeling, and relativistic Fe-K lines, 
is beyond the scope of this paper, and will be reported else- 
where. 



The best-fit parameters of these models are listed in table 3. 
We note that by replacing the diskbb model with the compPS 
model, the best-fit value of Tin increased by 13%. This is be- 
cause a spectrum hardened by Comptonization was assumed in 
the compPS model, an intrinsic temperature Tin became lower, 
leading to an increase in r-i^. While the difference should be 
regarded as a systematic modeling uncertainty, we confirm that 
it does not affect our conclusion. The photon index of the hard 
tail (F = 2.13 ± 0.01) is smaller than the result obtained from 
the same Suzaku data by Reis et al. (2010) (F = 2.54^°;°?), 
who did not include the 15-20 keV and 45-50 keV bands in 
fitting the HXD/PIN spectrum. 




2 5 10 20 50 

Energy (keV) 



Fig. 5. (a) Suzaku response-unfolded uFu spec- 
trum of XTE J 1752-223 from the best-fit for 
wabs*(compPS+pexriv+gaus), and the residuals between the 
data and best-fit models for (b) wabs*(diskbb+powerlaw), 
(c) wabs*smedge*(diskbb+powerlaw+gaus), (e) 

wabs*(compPS+pexriv+gaiis). 

5. Discussion 

We have analyzed the MAXI/GSC, the Swift/XRT, and 
Suzaku data of XTE J 1752-223 in the high/soft state from 
MJD 55218.62 to 55281.63. The overall continuum spec- 
tra were reproduced by the MCD plus power-law model with 
an interstellar absorption. The innermost temperature Tin de- 
creased from ^0.7 keV to ^0.4 keV, while the innermost radius 
Tin remained constant. By assuming D = 3.5 kpc and i — 0°, 
the values of Tin derived from the MAXI/GSC, Swift/XRT, and 
Suzaku data are 41.0±2.1 km, 41.4±0.3 km, and42.9±0.9, re- 
spectively. The results from the three instruments are mutually 
consistent with each other within the statistical errors. For the 
following discussion, we employ 41 .0±2. 1 km for ri„ as deter- 
mined from the MAXI/GSC, although we also discuss the case 
when the Suzaku result with the compPS model (rin=47.2±0.8 
km) is adopted to take into account the possible systematic un- 
certainty. 
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Since the value of rj„ thus estimated is significantly larger 
than those found in luminous low-mass X-ray binaries (10 
km; Mitsuda et al. 1984), the black hole interpretation of 
XTE J 1752-223 (Paper I; Munoz-Darias et al. 2010) is con- 
siderably reinforced. Furthermore, its constancy allows us to 
identify it with the ISCO in the high/soft state. The stable rjn 
at 41 .0 km is supposed to reach the ISCO in the high/soft state. 
We note that this rjn is an "apparent" innermost radius, and the 
"realistic" innermost radius (i?in) should be estimated as i?in= 
^K^Hn where the spectral hardening factor/* is 1.7 (Shimura & 
Takahara 1995) and correction factor for the boundary condi- 
tion ^ is 0.412 (Kubota et al. 1998). The value of k has been 
confirmed in recent work on disk models (e.g., Kubota et al. 
2010; Done & Davis et al. 2008). When the central object is 
assumed to be a non-spinning black hole (i.e. Schwarzschild 
black hole), R-m should equal to 6i?g. Then, the black-hole 
mass is estimated as 

M=^.5.51±0.28(^)(cos.)-M,. (1) 

Using the inclination angle i <49° obtained from radio obser- 
vations (Miller- Jones et al. 201 1), the black hole mass and dis- 
tance are constrained in the shadowed area shown in figure 6. If 
we assume the distance to be 3.5 kpc, the mass will be 5.2-7. 1 
Mq, and if i:>=10 kpc, 15.0-20.4 Mq. 

To validate our method to estimate the black hole mass by 
employing the MCD model which is a simplified approxima- 
tion of a true disk spectrum, we also performed a spectral 
fit to the Suzaku spectra using the kerrbb model (Li et al. 
2005). This model calculates the X-ray spectrum of a rel- 
ativistic accretion disk around a rotating black hole by tak- 
ing into account the innermost boundary condition and the 
effects of self-irradiation of the disk. Here we fixed the dis- 
tance at _D=3.5 kpc, spin parameter at a=0, and color correc- 
tion factor at k=\.1 for consistency with the previous discus- 
sion. The fit with the wabs*(kerrbb+pexriv) model was not ac- 
ceptable ix^/dof = 215.5/135), leaving similar residuals seen 
in figure5(b). Thus, we employed the simpl model (Steiner et 
al. 2009) with a steep photon index fixed at 7.4 to approx- 
imately represent weak Comptonization of the disk emission; 
the model became wabs*(simpl*kerrbb+pexri\'+gaus), which 
was found to give an acceptable fit (x^ /dof = 129.5/134). The 
black hole mass derived for i=0° and i =49° is 4.98t°;25 
and 8.43to:42 M©, respectively, which differs by --20% from 
the mass presented in equation (1). 

From past observations of BHCs and neutron stars 
Maccarone (2003) pointed out that the state transition from 
the high/soft state to the low/hard state occur at 1-4% (cen- 
tered at 2%) of the Eddington luminosity ^. Judging from the 
change of the photon index, XTE J 1752-223 went back to the 
low/hard state around MJD 55292 (DatalD 49). The bolo- 
metric flux at the transition to be Fti.ans=(6.68±0.83)x 10^^ 
erg cm~^ s~^ assuming a cutoff power-law is continuum with 
r = 1.8 and a cutoff energy of 200 keV. Using the rela- 
tion 0.01 < £trans/iodd <0.04 and Ltrans = 4:TrD'^Fbo\, we 
can further constrain the distance and black hole mass of 
XTE J 1752-223 as illustrated with the hatched region in fig- 

^ Lcdd = l-5x 10^** -p^erg s^^ for the solar abundances 




Distance (kpc) 



Fig. 6. Observational constraints on the mass and distance of 
XTE J 1752-223. The shadowed region is specified by the innermost 
radius derived from the MAXI/GSC spectra in combination with the 
inclination angle «<49° from the radio observation (Miller- Jones et 
al. 2011). Each solid line with two dashed lines indicates the best fit 
parameter and 90% confidence ranges of flux and rjn. The hatched 
region is derived assuming the empirical relation that the soft-to-hard 
transition occurs at 1~4 %L^cid- 

ure 6. Thus we can obtain D=2.5-\% kpc and A/=3.7-36 Mq. 
When we instead use the Suzaku results, L)=3.0-20 kpc and 
A/=5.3-44 Mq (based on the wabs*(compPS+pexriv+ gaus) 
model) or D=2.3-22 kpc and A/=3.1-55 Mq (based on the 
wabs*(simpl*kerrbb+pexriv+gaus)) are derived. 

Since we obtain M >3Mq, our method of estimating M 
assuming a black hole is self consistent. The conclusion that 
the compact object in XTE J1752-223 is a black hole is robust 
against the assumptions that the rotational axis of the accre- 
tion disk and jet axis is exactly aligned and that the black hole 
is non-spinning, since the extreme conditions i = 0° and a=0 
yield the smallest mass estimate. The values of M=9.6±0.9 
A/0 and _D=3.5±0.4 kpc estimated by Shaposhnikov et al. 
(2010) are outside the hatched region. If the black hole in XTE 
J 1752-223 is spinning and/or the disk has a larger inclination 
than 49°, then the discrepancy could be solved. 

6. Conclusions 

Using data obtained by the MAXI/GSC, the Swift/XRT, 
and Suzaku, we have performed X-ray spectral analysis of the 
black hole candidate XTE J 1752-223 in the high/soft state. As 
commonly seen in BHCs, the innermost radius remained con- 
stant in this state during the continuous observation with the 
MAXI/GSC. The results using the MCD plus power-law model 
were consistent between the three observatories. Detailed 
spectral modeling using Suzaku data suggests a possibility 
that the MCD emission may be slightly Comptonized, which 
could explain apparently broad iron-K features. Assuming 
that the innermost radius of the disk corresponds to the ISCO 
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and the black hole is non-spinning, we estimated the mass 
of the central object as 5.51±0.28 Mq (3-^^) (cosi)"i 
with corrections for the stress-free inner boundary condition 
and color hardening factor of 1.7. Furthermore, using the 
observational results that the inclination angle ; is smaller than 
49° from radio observations and the "soft-to-hard" transition 
occurs at 1-4 % Eddington luminosity, the distance and black 
hole mass of XTEJ1752-223 were constrained as 3.1-55 
Mq and 2.3-22 kpc, respectively, based on an analysis of the 
Suzaku spectra with a relativistic accretion disk model. Thus 
we can conclude that XTE J 1752-223 is likely a black hole. 

This research has made use of MAXI data provided by 
RIKEN, JAXA and the MAXI team. Suzaku data is pro- 
vided by a collaborative mission between the space agencies 
of Japan (JAXA) and the USA (NASA). We also thank the 
Swift team for their observation. This work made use of 
data supplied by the UK Swift Science Data Centre at the 
University of Leicester. This research was partially supported 
by the Ministry of Education, Culture, Sports, Science and 
Technology (MEXT), Grant-in-Aid No. 19047001, 20041008, 
20540230, 20244015 , 20540237, 21340043, 21740140, 
22740120. One of the authors (S. N.) is grateful to a grant 
from the Hayakawa Satio Fund awarded by the Astronomical 
Society of Japan. 
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Table 1. Best fit parameters of MAXI observations. 
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1 . 1 11^ J j^ 




55224, 


.84- 


-55225, 


.Oj 


7 07 


z.zU(,nxj 


A-l-O C 
Z.OitU.O 


O A^-i-O O/l 
U.0j±U.U4 


70 c+B.e 
32.5_5 2 


8. 


.0±0.4 


1 1 ^/'7^^ 
1 . 1 j^jj ) 


Zi 


55226, 


.30- 


-55227, 


.i i 


J. J 1 


z.zUt^nxj 


7 A^O 7 


A Al -i-O O/l 
U.01±U.U4 


7A Q+7-7 


7, 


.4±0.4 


A7/'7Q^ 

u.Oj^jy ) 


LL 


55227, 


.37- 


-55228, 


7Q 

iy 


7 ^7 
J.J J 




1 A ^O 7 
1 .4±U. / 


O AO-i-O 07 
U.OZitU.Uj 


78 0+6-2 
jO.U_g 2 


8. 


.2±0.4 


1 1 7/'77^ 

1 . 1 j^j J ) 


ZJ 


55228, 


.45- 


-55229, 


1A 
. j4 


J. J J 


z.zU(,nxj 


1 ^^O 7 
1 . J±U. / 


O AO-i-O O/I 
U.0U±U.U4 


78 O+'^'-O 
3S.2_5 2 


7, 


.1±0.4 


1 ^Q/'7o^ 
1 . ly^juj 


O /I 

z4 


55229, 


.41- 


-55230, 


.36 


1 oo 
3.0Z 


z.zU(nx) 


1 O ^A A 
1.0±0.6 


A C C _l_ A A7 
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CI -7 + 10.5 
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7. 
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1 0C/70\ 
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o^ 

ZJ 
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.55- 
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o 7^o n 

L. I±\J. 1 


O <7-i-A O/l 
U.J /±U.U4 


A A T+10.1 

44.3_8o 


6, 


.7±0.4 


1 OQ/'7/l^ 

i.zy(j4j 


o^; 
ZO 


55231 


.63- 


-55232, 


c\c\ 
.90 


A ir\ 
4.30 


z.zU(nx) 


Q+0-6 
^■^-0.7 


O ^A+0-04 


AO A + 10.1 

48.0_8.i 


6. 


.8±0.3 


A A7 /7 0\ 

O.y3(3o) 


OT 

Z / 


55237, 


.15- 


-55237, 


AO 
.OZ 


O AA 

z.oo 


z.zU(ilxj 


1 <^A T 


A AA-i-A A^ 

U.oU±U.Uj 


7A A+9-3 


6 


.0±0.4 


A 70/'0A\ 

u. /y(zo) 


o o 
ZO 


55237, 


.87- 


-55238, 


A C 

.45 


O O /I 

z.o4 


z.zU(nx) 


1 *7^A A 
1. /±0.6 


A C O _l_ A AC 

0.5o±0.05 


7Q 7 + 10.8 

39.3„8 1 


5, 


.8±0.4 


A AA/0-7\ 

o.y6(z/) 


oo 

zy 


55238, 


.82- 
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2, 


.3±0.1 
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44 
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.28- 
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2.20(fix) 


0.5±0.2 


0.46±0.03 
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-'l^--'-10.5 


2, 


.0±0.1 
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45 


55277, 


.23- 


-55281, 


.63 


3.98 


2.20(fix) 


0.8±0.5 


47+0-09 

^■^'-0.08 


39 0+3^-2 

--'^■^-16.6 


1, 


.4±0.2 


0.94(14) 


46 


55285, 


.44- 


-55286, 


.97 


2.34 


2.20(fix) 


^■1-0.5 


31+0-20 

"■--'1-0.14 


118.7l?i™7^ 





.5±0.2 


0.46(9) 


47 


55287, 


.41- 


-55288, 


.05 


2.42 


Z.OJ_Q 4Q 


2.0±0.4 










0.39(9) 


48 


55288, 


.43- 


-55291, 


.16 


8.98 


2 19+0-21 


2.0±0.2 










0.58(32) 


49 


55291 


.35- 


-55292. 


.49 


3.41 




o 7+0.5 
^- '-0.4 










0.61(14) 



Note^ 

* In unit of 10"'' erg s"! cm-^ (2-20 keV). 
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Table 2. Best-fit parameters of Swift/XRT observations. 



ObsID 



MJD 

start 



Exposure 
(sec) 



Nh 

(xl022 cm-2) 



power-law 
flux* 



Tin 

(keV) 



(km) 



Disk 
flux* 



38.4fTy 

38.0±1.0 

44 3+0.8 

^■•'-0.7 

42.2±1.1 
40.9^1^ 
41.7±l'.l 
45.3±1.4 
41.7±1.4 
41.9±0.7 
43.4±1.1 

+ 1.6 
-1.5 

36.9ll:J 

-'^■■'-0.9 

43.3±1.3 

+ 1.2 
-1.1 
+3.1 
-2.8 



Xi /dof 



00031532010 
00031532011 
00031532012 
00031532014 
00031532015 
00031532016 
00031532018 
00031532019 
00031532020 
00031532021 
00031640001 
00031640002 
00031640003 
00031640004 
00031640005 
00031640006 
00031640007 



55231.6 
55233.7 
55238.9 
55242.0 
55244.0 
55246.2 
55248.0 
55250.1 
55251.3 
55252.6 
55257.6 
55260.7 
55266.8 
55269.8 
55271.8 
55276.9 
55279.8 



1374.3 

1595.0 

2360.7 

1421.4 

1123.8 

1411.8 

835.6 

1144.3 

3026.5 

1524.0 

890.9 

909.9 

1209.8 

1184.1 

1214.6 

977.8 

128.3 



Notes. 

Photon indices was fixed to 2.20. 

• In unit of 10"^ erg s"! cm-^ (2-20 keV). 



0.56±0.01 

0.60±0.01 
0.62±0.01 
0.56±0.01 
0.55±0.01 
0.59±0.01 
0.63±0.01 
0.61±0.01 
0.61±0.01 
0.57±0.01 
0.60±0.01 
0.61±0.01 
0.61±0.01 
0.67±0.01 
0.66±0.01 
0.57±0.01 

"■-'^-0.02 



2.34±0.23 
2.40±0.19 
1.26±0.11 
2.41±0.19 
1.07±0.17 
2.80±0.18 
1.29±0.16 
0.47±0.12 
0.87±0.07 
0.65±0.09 
0.69±0.11 
0.25±0.05 
0.17±0.06 
0.32±0.04 
0.33±0.04 
0.34±0.04 
±0.12 



0.60±0.01 
0.59±0.01 
0.58±0.01 
0.57±0.01 
0.56±0.01 
0.57±0.01 
0.55±0.01 
0.54±0.01 
0.55±0.01 
0.53±0.01 
0.52±0.01 
0.52±0.01 
0.50±0.01 
0.50±0.01 
0.48±0.01 
0.48±0.01 
0.47±0.01 



42.5 



41.5 
35.8 



5.66±0.20 
5.39±0.14 
6.39±0.11 
5.21±0.15 

4 Q1+0.17 

4 7Q+0.13 

5 01+" " 

•'■^"-0.14 

4.07±0.07 
3.80±0.06 
3.05±0.08 

9 9^+0.09 
^■^-'-0.08 

2.38±0.09 
2.00±0.05 
1.75±0.06 
1.67±0.04 

QQQ + O.ll 



1.68(224) 
1.29(238) 
1.50(244) 
1.64(234) 
1.22(191) 
1.56(241) 
0.95(202) 
1.22(177) 
1.62(238) 
1.50(193) 
1.22(174) 
1.34(180) 
1.25(161) 
1.61(197) 
1.41(181) 
1.34(177) 
1.00(108) 



Table 3. Best-fit parameters for tlie Suzaku observation. 



Model 


wabs diskbb 1 compPS^ 


powerlaw 1 pexriv^^ 


gaus** 


smedge^^ 




Tin (keV) 
rJn (km) 

r§ 


r 

norm* 


(keV) 
EW (eV) 


max 


1 d.o.f. 


wabs*( diskbb+powerlaw) 


0.62±0.02 0.53±0.01 
42.9±0.9 


2.29±0.06 
0.16±0.03 






230.1/136 


wabs *smedge *( diskbb+powerlaw+gaus ) 


0.65±0.02 0.526±0.003 
44.0±1.0 


2.22±0.05 
0.16±0.02 


"■-'^-0.07 
(ro+42 
"■'-41 




141.80/133 


wabs *( compPS+pexriv+gaus ) 


0.67±0.02 0.507l^;j;;jg 
47.2±0.8 
0.41±0.05 


2.13±0.01 
05+° °^ 


6 60+""^ 

"■""-0.07 

98±43 




131.4/133 



* In unit of 10^^ cm^-^. 

t To and covf^j^;. are fixed at 1 keV and at 1 , respectively. No reflection from the coinpPS 
component itself is considered. 
^D=3.5 kpc and j = 0° are assumed. 
§ Optical depth of the corona. 

II Fixed atQ = 2it, i = 25° and = 1000 erg cm/s. No high energy cutoff is set in the 
incident power law continuum. 

* A normalization in unit of photons cm^^ s^^ at 1 keV. 

* 1 (T line width is fixed at 10 eV. 

tt The edge energy and width, photo-electric cross-section are fixed at 7.11 keV, 10 keV 
and -2.67, respectively. 

it The maximum absorption factor at threshold. 
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Appendix 1. MAXI/GSC spectral calibration using the 
Crab Nebula 

We here show the status of the GSC energy response calibra- 
tion with Crab-nebula data to confirm the validity of the spec- 
tral analysis. We used event data with the process version 0.3b, 
which is screened with the operation high voltage of a nominal 
1650 V and the anode number of #0, #3, #4, and #5 whose posi- 
tional response had been well established. The event selection 
is same as that employed in the first performance verification 
in Sugizaki et al. (201 1). 

We here screened the event data with more severe conditions 
to verify the response calibration with a better accuracy. We 
selected events taken only during such a good scan transit that 
the source incident angle (pcoi is < 36° and the area for both the 
source and the background are completely covered. 

The calibration of the energy response was performed using 
the screened background and Crab Nebula data in the following 
steps; The energy-PHA relations, which depend on the detector 
position along the anode wires, are corrected by a gain factor 
for each output amplifier. The parameter was calibrated us- 
ing fluorescent lines in the background spectrum from Ti (4.5 1 
keV) and Cu (8.04 keV) as well as the calibration source ^^Fe 
(5.895 keV). We verified that the gain factor obtained with a 
0.1 % precision successfully reproduce gain position depen- 
dence within the discrepancy of 1% in RMS throughout the 
whole detector area. The threshold and the resolution of the 
Lower-Discriminator (LD) of each amplifier were then cali- 
brated using the Crab Nebula spectra. The LD parameters 
were tuned so that the results of the spectral fit to an absorbed 
power-law model agrees with those of the canonical values, 
photon index F = 2.1 and hydrogen absorption column density 
Nu ^ 0.35 X 10^2 cm-2 (i.e. Toor & Seward 1974, Kirsch et 
al. 2005). 

Figure A-1 shows the best-fit parameters to GSC Crab 
Nebula spectra. The data of multiple scan transits whose total 
exposure of >3.2 cm^ ksec were grouped so that each spec- 
trum has enough photon statistics to constrain parameters of 
a power-law model. Table A-1 summarizes the results. The 
derived parameters agree with the canonical values within the 
statistical errors. 



Table A-1. Best-fit values and variations of absorbed power-law model 
to GSC Crab Nebula spectra 



Parameter 


Canonical 


Best fit 


rms* 


iVn (xlO^^ cm"^) 


0.35 


0.39 


0.25 


F 


2.10 


2.11 


0.05 


Ft 




2.11 


0.03 


FluX2-20koV t 


3.0±10% 


3.08 


0.06 



Notes. 

* Root mean square of the best-fit value 
t In unit of 10^** erg s^^cm^^ 

* A^H is fixed at 0.35x10^2 cm-^ 




55150 55200 55250 55300 55350 55400 



MJD 

Fig. A-1. Best-fit values and statistical errors of absorbed power-law 
model to GSC Crab Nebula spectra against the observation time; (top) 
hydrogen column density, (middle) photon index, (bottom) 2-20 keV 
absorbed flux. 



